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Abstract 

Gas  Turbine  Technologies  (GTT)  and  Politecnico  di  Torino,  both  located  in  Torino  (Italy),  have  been  involved  in  the  design  and  installation 
of  a  SOFC  laboratory  in  order  to  analyse  the  operation,  in  cogenerative  configuration,  of  the  CHP  100  kWe  SOFC  Field  Unit,  built  by  Siemens- 
Westinghouse  Power  Corporation  (SWPC),  which  is  at  present  (May  2005)  starting  its  operation  and  which  will  supply  electric  and  thermal 
power  to  the  GTT  factory.  In  order  to  take  the  better  advantage  from  the  analysis  of  the  on-site  operation,  and  especially  to  correctly  design 
the  scheduled  experimental  tests  on  the  system,  we  developed  a  mathematical  model  and  run  a  simulated  experimental  campaign,  applying  a 
rigorous  statistical  approach  to  the  analysis  of  the  results. 

The  aim  of  this  work  is  the  computer  experimental  analysis,  through  a  statistical  methodology  (2k  factorial  experiments),  of  the  CHP  100  per¬ 
formance.  First,  the  mathematical  model  has  been  calibrated  with  the  results  acquired  during  the  first  CHP  100  demonstration  at  EDB/ELSAM 
in  Westerwoort.  After,  the  simulated  tests  have  been  performed  in  the  form  of  computer  experimental  session,  and  the  measurement  uncertain¬ 
ties  have  been  simulated  with  perturbation  imposed  to  the  model  independent  variables.  The  statistical  methodology  used  for  the  computer 
experimental  analysis  is  the  factorial  design  (Yates’  Technique):  using  the  ANOVA  technique  the  effect  of  the  main  independent  variables 
(air  utilization  factor  C70X,  fuel  utilization  factor  t/F,  internal  fuel  and  air  preheating  and  anodic  recycling  flow  rate)  has  been  investigated  in 
a  rigorous  manner.  Analysis  accounts  for  the  effects  of  parameters  on  stack  electric  power,  thermal  recovered  power,  single  cell  voltage,  cell 
operative  temperature,  consumed  fuel  flow  and  steam  to  carbon  ratio.  Each  main  effect  and  interaction  effect  of  parameters  is  shown  with 
particular  attention  on  generated  electric  power  and  stack  heat  recovered. 

©  2005  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

The  EOS  Project  -  partners:  Politecnico  di  Torino,  Gas 
Turbine  Technologies  (GTT,  Siemens  group),  Hysylab 
(Hydrogen  System  Laboratory)  of  Environment  Park  and 
Regione  Piemonte  -  aim  to  create  the  main  node  of  a  regional 
fuel  cell  generator  network.  As  a  first  step,  the  Pennsylvania- 
based  Stationary  Fuel  Cells  division  of  Siemens- 
Westinghouse  Power  Corporation  (SWPC)  supplied  GTT 
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with  a  CHP  100  kWe  SOFC  Field  Unit,  fuelled  by  natural 
gas  with  internal  reforming.  The  fuel  cell  is  connected  to  the 
electricity  national  grid  and  provides  part  of  the  industrial 
district  energy  requirement.  We  referred  to  some  papers  to 
acquire  some  information  before  designing  the  cogeneration 
plant  of  the  EOS-CHP  100  and  we  took  advantage  on  this  in 
planning  the  EOS-CHP  100  experimental  campaign. 

Most  publications  about  SOFC  mainly  focus  on  single 
cell  or  involve  aspects  of  the  balance  of  plant  (BoP)  of  a 
SOFC  plant.  These  papers  discuss  fundamental  scientific 
topics  of  the  main  cell  elements  and  of  the  BoP  devices: 
electrochemical  reactions,  fuel  processing  reactions,  cell 
component  overpotentials,  materials  and  degradation  of 
performances.  Most  papers  describe  the  development  of 
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Nomenclature 

A 

symbol  of  the  factor  Uox 

B 

symbol  of  the  factor  I/p 

C 

symbol  of  the  factor  7p 

CHP 

Cogeneration  heat  and  power 

D 

symbol  of  the  factor  Tox 

E 

symbol  of  the  factor  / 

Ea 

activation  energy  of  anode  (kJ  mol-1) 

Ec 

activation  energy  of  cathode  (kJ  mol-1) 

FVev 

open  circuit  voltage  (V) 

/ 

/-value  of  the  ANOVA  analysis 

F 

Faraday  number  (C  mol-1) 

G 

mass  flow  (kg  s-1) 

Gf 

fuel  mass  flow  (kg  s-1) 

h 

specific  enthalpy  (kJkg-1) 

h 

molar  enthalpy  (kJ  mol-1) 

i 

current  density  (A  cm-2) 

fi,a/c 

limiting  current  density  at  anode  and  cathode 
(A  cm-2) 

k 

number  of  factors  in  the  factorial  analysis 

Ka 

pre-exponential  coefficient  for  anode  in 
activation  overpotential  (Am-2) 

Kc 

pre-exponential  coefficient  for  cathode  in 
activation  overpotential  (Am-2) 

Kp 

equilibrium  constant  of  a  chemical  reaction 

Kpreforming  equilibrium  constant  of  the  reforming 
reaction 

KPshift 

equilibrium  constant  of  the  water  shift  reaction 

ma 

anode  activation  overpotential  exponential 

mc 

cathode  activation  overpotential  exponential 

n 

number  of  electrons  involved  in  an  electrode 

reaction 

nc 

number  of  tubular  cells  in  series  in  the  stack 

^ch4, inlet  number  of  methane  moles  at  the  inlet  of  the 
fuel  mixer  (mol  s-1) 

^H20,rec 

number  of  steam  water  moles  recirculated 
from  anode  to  fuel  mixer  (mol  s-1) 

n\ 

total  number  of  moles  at  fuel  mixer  inlet  flow 
(mols-1) 

no2,  consumed  number  of  oxygen  moles  consumed  inside 
the  stack  (mol  s-1) 

n02,  inlet 

number  of  oxygen  moles  at  the  inlet  of  the 
stack  (mol  s-1) 

Nk 

molar  chemical  reaction  rate  (mol  s-1) 

P 

pressure  (bar) 

Pstack 

stack  operating  pressure  (bar) 

R 

universal  gas  constant  (J  mol-1  K-1) 

Ra 

resistance  equivalent  of  the  anode  activation 
overpotential  (£2) 

Rc 

resistance  equivalent  of  the  cathode  activation 
overpotential  (£2) 

r 

resistance  of  a  layer  in  ohmic  overpotential 
(£2  cm2) 

S 

single  cell  active  area  (cm2) 

S/C 

steam-to-carbon  ratio 

SOFC 

solid  oxide  fuel  cell 

T 

temperature  (K) 

Tf 

fuel  pre-heating  temperature  (°C) 

Tox 

air  pre-heating  temperature  (°C) 

/stack 

stack  operating  temperature  (°C) 

t 

thickness  of  a  layer  in  ohmic  overpotential 
(cm) 

Uf 

fuel  utilization  factor 

Gox 

air  utilization  factor 

Vc 

single  cell  voltage  (V) 

w 

effect  of  a  factor  evaluated  through  the  Yates’ 
Technique 

^el 

stack  electrical  power  (kW) 

Wel,c 

single  cell  electrical  power  (W) 

V 

moles  of  CH4  which  react  (mol  s-1) 

J 

moles  of  CO  which  react  (mol  s-1) 

z 

moles  of  H2  which  react  (mol  s-1) 

Greek  letters 

AG 

change  in  Gibbs  free  energy  in  a  reaction 
(J  mol-1) 

^ech 

heat  of  electrochemical  reaction  (kW) 

0p 

heat  for  fuel  pre-heating  (kW) 

^ox 

heat  for  oxidant  (air)  pre-heating  (kW) 

^  ref 

heat  of  steam  reforming  reaction  (kW) 

shift 

heat  of  water  shift  reaction  (kW) 

^th 

stack  recovered  heat  power  (kW) 

P 

fraction  of  anode  gas  recycling 

7?act,a/c 

activation  overpotential  at  anode  and  cathode 
(V) 

7?conc,a/c 

concentration  overpotential  at  anode  and 
cathode  (V) 

?7  ohm 

ohmic  overpotential  (V) 

Pa. 

resistivity  of  anode  layer  (£2  cm) 

Pc 

resistivity  of  cathode  layer  (£2  cm) 

Pc  1 

resistivity  of  electrolyte  layer  (£2  cm) 

Pint 

resistivity  of  interconnection  (£2  cm) 

modeling  studies  of  single  solid  oxide  fuel  cell  or  BoP,  or 
the  results  of  experimental  analysis  on  single  tubular  or 
planar  cell.  Some  other  papers  are  devoted  to  the  activity 
concerning  SOFC  power  system’s  demonstrations. 

In  literature,  many  papers  take  into  account  a  parametric 
analysis  of  a  SOFC  plant  by  changing  the  stack  operating 
conditions.  The  aim  of  these  approaches  is  to  evaluate  the 
performance  of  the  plant  by  changing  the  BoP  configura¬ 
tion  internal  or  external  fuel  processing,  combined  heat  and 
power  generation  plants  or  hybrid  SOFC/GT  plants  [1-10]) 
or  the  generator  operation  and  the  fuel  feeding.  The  effects 
of  the  modification  of  these  variables  were  not  analysed  in 
depth  with  statistical  methodologies,  reducing  to  a  qualitative 
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description  of  the  effects  on  the  performance  of  the  simulated 
stack. 

Some  fundamental  topics  about  SOFC  are  also  presented 
by  experimental  characterization.  Some  papers  point  out 
the  evaluation  of  the  electrochemical  behavior  of  SOFC 
anodes  [11-16].  Other  papers  show  experimental  analysis  in 
order  to  investigate  degradation  phenomena  at  cell’s  compo¬ 
nents  under  typical  operating  conditions;  both  macroscopic 
and  micro-structural  degradations  are  investigated  like  car¬ 
bon  deposition,  delamination,  nickel  agglomeration,  thermal 
cycling  [16-29] ;  some  studies  focus  on  mechanical  properties 
of  SOFC  materials  by  evaluating  creeps  or  curvature  under 
different  test  conditions  [29-31]. 

1.1.  Modeling  approach 

Modeling  of  SOFC  generally  involves  two  main 
approaches:  the  development  of  numerical  tools  with  the  aim 
to  deeply  describe  chemical,  electrochemical  and  thermal- 
fluid  dynamics  processes  taking  place  in  a  single  cell  [32-42] 
or  the  development  of  zero-dimensional  models  with  the  aim 
of  modeling  a  SOFC  generator  integrated  in  a  more  complex 
power  generation  plant.  This  approach  allows  to  perform  sim¬ 
ulated  generator  steady-state  operation  with  the  possibility  of 
taking  suggestions  about  optimal  BoP  configuration  or  gener¬ 
ator  operation.  In  the  development  of  our  model,  we  referred 
mostly  to  the  second  approach.  Some  of  the  main  considered 
papers  are  discussed. 

In  [1],  the  mathematical  simulation  of  a  planar  solid 
oxide  fuel  cell  (SOFC)  is  presented.  It  accounts  for  three- 
dimensional  and  time- dependent  effects.  Internal  methane- 
steam  reforming  and  recycling  of  the  anode  gas  are  also 
considered.  In  [2],  a  mathematical  model  of  a  tubular  solid 
oxide  fuel  cell  (SOFC)  is  presented.  The  complete  electro¬ 
chemical  and  thermal  factors  are  accounted  for  in  a  rigorous 
manner.  Comparison  is  made  with  single  cell  test  data  from 
Westinghouse.  In  [3],  the  aim  of  the  work  is  to  investigate 
of  internal  reforming  solid  oxide  fuel  cell  and  gas  turbine 
combined  cycles.  A  mathematical  model  has  been  developed 
that  simulates  fuel  cell  steady-state  operation,  the  model  has 
been  used  for  a  complete  parametric  analysis;  after  in  [4] 
both  an  exergy  and  a  thermoeconomic  analysis  of  the  pro¬ 
posed  cycles  have  been  carried  out;  recently  in  [5]  presented 
a  work  dealing  with  the  design  and  off-design  performance 
evaluation  of  an  anodic  recirculation  system  based  on  ejector 
technology  for  solid  oxide  fuel  cell  hybrid  applications.  In 
[6],  a  thermodynamic  model  of  a  tubular  SOFC  stack  is  pro¬ 
posed.  It  is  then  calibrated  on  the  available  data  for  a  recently 
demonstrated  tubular  SOFC  prototype  plant  (Westerwoort; 
CHP  100  Siemens).  In  [7],  a  detailed  BoP  has  been  investi¬ 
gated.  An  energy  and  economic  analysis  of  a  decentralized 
natural  gas-fuelled  SOFC-power  plant  in  the  range  of  200  kW 
capacity  is  carried  out.  All  the  calculations  start  from  a  basic 
plant  concept  with  a  simple  flowsheet  and  a  basic  param¬ 
eter  set  of  SOFC  operation  and  economic  data.  Changes  in 
costs  of  electricity  and  plant  efficiency  are  determined  for  the 


variation  of  the  cell  operation  parameters.  This  includes  the 
influence  of  the  air  temperature  increase  in  the  stack,  degree 
of  internal  reforming,  cell  voltage  and  fuel  utilization;  in  the 
Part  II  of  the  paper  [8],  the  flowsheet  of  the  SOFC  BoP  is 
changed  considering  gas  recycling  by  blowers  or  jet  boost¬ 
ers. 

Some  articles  presented  recently  SOFC  technology 
demonstrations  and  SOFC  state  of  the  art  and  future  perspec¬ 
tives  [43-48].  In  particular,  in  [43]  the  Siemens’s  Field  Unit 
demonstration  program  including  the  first  SOFC  CHP  100 
demonstration  at  EDB/ELSAM  in  Westerwoort  (Holland), 
the  Southern  California  Edison  SCE  220-kW  and  pressurized 
SOFC/gas  turbine  (PSOFC/GT)  power  system  are  described. 

1.2.  Aims  of  the  paper 

Some  authors  finalized  their  papers  in  a  parametric 
analysis  of  a  SOFC  system.  In  [3, 4, 6-8],  the  effects  of 
cell  temperature,  cell  pressure,  recirculation  ratio  of  anode 
exhaust,  utilization  ratios  of  fuel  and  air,  and  average  current 
density  on  both  the  electric  and  thermal  performances  of  a 
solid  oxide  fuel  cell  were  calculated  under  typical  operating 
conditions.  The  effects  of  these  operation  independent 
variables  were  not  analysed  in  depth  with  statistical  method¬ 
ologies,  reducing  to  a  description  of  a  parametric  analysis. 

Now,  we  are  planning  the  experimental  campaign  of  the 
EOS  Project  and,  in  order  to  take  the  better  advantage  from 
the  analysis  of  the  on-site  operation,  and  especially  to  cor¬ 
rectly  design  the  scheduled  experimental  tests  on  the  system, 
we  developed  a  model  of  the  EOS-CHP  100  and  run  a  simu¬ 
lated  experimental  campaign,  applying  a  rigorous  statistical 
approach  to  the  analysis  of  the  results. 

The  aim  of  our  paper  is  therefore  the  experimental  anal¬ 
ysis,  through  a  statistical  methodology,  of  the  single  and 
combined  effects  of  the  main  stack  operation  independent 
variables  on  the  cogenerative  (electricity  and  heat)  perfor¬ 
mance  of  a  tubular  SOFC  stack  integrated  system. 

The  tests  have  been  performed  in  the  form  of  a  computer 
experimental  session,  based  on  a  simulation  model  of  the 
stack.  The  measurement  uncertainties  are  simulated  with  per¬ 
turbation  imposed  to  some  independent  variables. 

The  statistical  methodology  used  for  the  computer  exper¬ 
imental  data  analysis  has  been  the  factorial  design  (Yates’ 
Technique):  using  analysis  of  variance  (ANOVA)  we  have 
analysed  the  significance  of  the  main  operation  independent 
variables  (factors:  air  utilization  factor,  fuel  utilization  fac¬ 
tor,  fuel  pre-heating  temperature,  air  pre-heating  temperature, 
fraction  of  anode  gas  recycling)  considering  their  single  and 
combined  effects  on  the  electric  and  thermal  power  which 
could  be  recovered  from  the  stack  system. 

2.  Description  of  SOFC  stack  and  BoP 

The  EOS-CHP  1 00  Field  Unit  is  the  first  field  unit  to  utilize 
the  commercial  prototype  of  air  electrode  supported  cell  and 
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Fig.  1.  Picture  of  the  EOS-CHP  100  SOFC  Field  Unit  installed  at  the  GTT 
factory. 

in- stack  reformers  (Fig.  1).  In  the  design  used  in  the  EOS- 
CHP  100  SOFC  Field  Unit,  the  cell  components  are  deposited 
in  the  form  of  thin  layers  on  the  air  electrode  (cathode)  tube 
(air  electrode- support  AES)  closed  at  one  end.  The  size  of  the 
complete  tube  is  150  cm  length  and  2.2  cm  diameter  with  an 
external  surface  (anode  side)  of  834  cm2  [43-47].  The  next 
level  of  fabrication  hierarchy  after  the  cell  is  the  cell  bundle, 
which  consists  of  a  24-cell  array  arranged  as  eight  cells  in 
electrical  series  by  three  cells  in  electrical  parallel.  Four  cell 
bundles  are  connected  in  series  to  form  a  bundle  row,  and 
12  bundle  rows  are  aligned  side  by  side,  interconnected  in 
serpentine  fashion  with  an  in- stack  reformer  between  each 


bundle  row  (4  bundles  in  12  rows,  for  a  total  of  1 152  single 
tubes).  An  advantage  of  all  cathode-supported  single  cells  is  a 
three-dimensional  penetration  structure  in-between  cathode 
and  electrolyte  that  is  achieved  during  electrolyte  deposition. 
This  leads  to  a  decrease  of  cathodic  polarization  losses  [45]. 

A  detailed  description  of  materials  and  geometry  charac¬ 
terizing  a  single  solid  oxide  fuel  cell  tube  could  be  found  in 
literature  [43-47]. 

Considering  the  total  volume  of  the  EOS-CHP  100  SOFC 
Field  Unit,  just  the  25%  is  occupied  by  the  stack.  The  rest 
of  the  system  is  composed  by  the  balance  of  plant  (BOP), 
including  five  major  skids:  Generator  Module,  Electrical 
Control  System,  Fuel  Supply  System  (FSS),  Thermal  Man¬ 
agement  System  (TMS)  and  Heat  Export  System  (HES).  The 
simplified  flow  schematic  of  the  BOP  of  the  EOS-CHP  100 
SOFC  Field  Unit  is  shown  in  Fig.  2. 

The  stack  is  integrated  with  a  catalytic  pre-reformer  for  the 
natural  gas-fuel,  and  a  combustor  where  the  post-combustion 
of  the  depleted  fuel  and  air  exiting  the  cell  occurs.  Before 
entering  the  stack,  the  fuel  flow  is  desulfurized  and  then 
sent  to  the  ejectors  in  the  stack.  Part  of  the  anode  exhaust 
gases  is  recirculated  in  the  fuel  stream  (to  supply  water  for 
the  fuel  internal  reforming  reaction).  Since  the  generator 
operates  at  high  temperature  (around  970  °C),  a  heat  flow 
from  the  stack  exhausts  can  be  recovered  at  the  outlet  of  the 
plant.  Considering  the  assessment  of  the  BoP  of  the  EOS- 
CHP  100,  the  exhausts  leave  the  generator  skid  at  around 
850-900  °C.  After,  in  the  thermal  management  system  (TMS) 
they  exchange  energy  with  the  incoming  process  air  (in  the 
low  temperature  recuperator  (LTR)  and  in  the  high  temper¬ 
ature  recuperator  (HTR));  at  the  TMS  outlet,  at  nominal 
operating  condition  (generator  current  500  A  and  Up  =  0.85) 


Fig.  2.  Simplified  flow  schematic  of  the  balance  of  plant  of  the  EOS-CHP  100  SOFC  Field  Unit. 
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the  exhausts  temperature  is  around  250  °C.  In  the  heat  export 
system,  the  heat  is  recovered  from  the  stack  exhausts  in  order 
to  heat  a  pressurized  water  flow  which  will  be  used  in  the 
heating  system  of  the  GTT  factory  (during  Winter)  or  to  pro¬ 
duce  cooling  power  in  an  absorption  chiller  (during  Summer). 
Therefore,  the  present  configuration  of  the  EOS-CHP  plant  is 
a  trigeneration  energy  system.  The  auxiliaries  are:  process  air 
blowers,  air  filters,  purging  gas  vessels,  start-up  and  control 
systems  [43-45]. 

In  the  nominal  operation  point  (7=500  A,  corresponding 
to  0.2  A  cm-2),  the  DC  electric  efficiency  is  53.25%,  while 
the  AC  electric  efficiency  is  46%,  while  the  cogenerative 
efficiency  (AC  electric  +  recovered  thermal  energy)  is  in  the 
order  of  74-75%. 


3.  SOFC  stack  model 

Modeling  the  SOFC  stack  involves  the  electrochemical 
and  chemical  models.  The  electrochemical  model  refers  to 
the  polarization  curve.  The  chemical  model  considers  the 
reactions  taking  place  in  the  stack:  fuel  processing  reactions 
(steam  reforming,  water  shift). 

The  subsequent  balance  of  plant  will  integrate  the  stack 
model  with  a  thermal  model  (taking  into  account  the  energy 
balance  equations)  and  a  chemical  species  mass  balance  for 
the  balance  of  plant. 


3.1.  Electrochemical  model 


Table  1 

Equations  of  the  electrochemical  model 
Activation  overpotential  rjac t>a/c  [1] 

(\  tnc 

fr)  exP(-Jr) 

(3)Anode  V  =  .Ka  •  ■  (^-)  exp(-|f=) 

Ohmic  overpotential  r\ 0hm  [2,32] 

(4)^ohm  =  r  j  =  pjtj 

i 

(5)Anodepa  =  0.008114  exp 

(6)Electrolyte  Pel  =  0.00294  exp 

(7)Cathodepc  =  0.00298  exp 
^Interconnection  p-m  =  0.1256  exp 

Concentration  overpotential  /ponc.a/c 

(9)Anode/ cathode  rjc onc,a/c  =  fp  In  (  1  -  r1-  ) 


Constants 

Activation  energy  of  anode  1  lOkJ  mol-1 
Activation  energy  of  cathode  160kJ  mol-1 
Pre-exponential  coefficient  for  anode  Ka  =  2.13  x  108Am  2 
Pre-exponential  coefficient  for  cathode  Kc  =  1 .49  x  1010  Am 
Anode  activation  overpotential  exponential  ma  =  0.25 
Cathode  activation  overpotential  exponential  mc  =  0.25 


These  reactions  are  considered  at  equilibrium.  For  both  reac¬ 
tions,  the  equilibrium  constant  is  calculated  as: 


In  Kp(T) 


-AG(T,  p°) 
RT 


(12) 


The  polarization  curve  for  a  solid  oxide  fuel  cell  is 
expressed  by  the  Eq.  (1): 

Fc  —  Frev  —  ?7act,a/c  —  Oohm  —  7?conc,a/c  (1) 

where  Vc  is  the  single  cell  voltage,  Frev  the  open  circuit 
voltage,  ?7act,a/c  the  activation  overpotential,  r/ohm  the  ohmic 
overpotential  and  0conc,a/c  is  the  concentration  overpotential. 

Most  of  the  considered  equations  are  suggested  by  current 
literature,  and  they  are  shown  in  Table  1.  In  particular,  the 
equations  modeling  the  resistivities  have  a  good  agreement 
with  the  available  experimental  data  in  the  range  of  current 
density  up  to  0.3  A  cm-2;  after  this  value,  the  models  are 
less  performing.  Up  to  now,  the  available  data  describing  the 
real  operation  of  the  stack  never  reached  values  higher  than 
0.3  A  cm-2. 

3.2.  Chemical  model 

The  chemical  model  accounts  for  internal  fuel  processing 
(steam  reforming,  water  shift)  reactions: 

Steam  reforming  :  CH4  +  H2O  3H2  +  CO  (10) 

(11) 


Eq.  (12)  can  be  used  integrated  by  a  correlation  linking 
directly  to  the  temperature  at  which  reaction  takes  place: 

log  Kp(T)  =  AT4  +  BT3  +  CT2  +  DT  +  E  (13) 

The  parameters  A-E  of  the  JANAF  tables  have  been  found 
in  literature  [3]. 

4.  BoP  model 

The  BOP  model  accounts  for  devices,  which  define  the 
SOFC’s  plant.  These  components  are:  blowers,  desulfuriz- 
ers,  purging  gas  tanks,  pre-reforming  reactor,  reformer,  heat 
exchangers  used  to  pre-heat  inlet  gases  (fuel  gas,  air),  heat 
exchanger  to  export  heat  by  wasted  stack  gases  (HES,  heat 
export  system)  and  post-combustion  chamber. 

The  hypothesis  for  the  BoP  model  are: 

•  steady- state  condition; 

•  zero-dimentional  aproach; 

•  cathodic  flow  composition:  mixed  flow  of  N2  and  O2; 

•  anodic  flow  composition:  mixed  flow  of  CH4,  CO2,  CO, 
H2,H20,  N2; 

•  internal  fuel  processing:  reforming  and  CO-shift  reactions 
at  equilibrium; 


Water  shift :  CO  +  H20  H2  +  C02 
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Table  2 

Equations  of  the  BoP  model 

Mixer  mass  balance  =  G3  +  G5  (15) 

Mixer  thermal  balance  G^h^  =  G3I13  +  G5I15  (16) 

Steam  to  carbon  ratio  (S/C)  S/C  =  /iH2°’rec  (17) 

7  «CH4, inlet 

Oxidant  consumption  n  o2inlet  =  n  o2consumed  /  %  (18) 

Heat  of  steam  reforming  reaction  0ref  =  x  •  (3/?h2  +  ^co  — 
hn2o  ~  ^ch4)  (19) 

Heat  of  water  shift  reaction  0shift  =  y  •  (^h9  +  ^co2  —  ^co  — 
^h2o)  (20) 

Heat  of  electrochemical  reaction  :  0ech  =  nc  •  I  • 

f  ~A/i2(^ack)  _  yc^  (21) 

Cell  electrical  power  :  VPei,c  =  J’cAVc  (22) 

Fuel  pre-heating  :  0F  =  Gp  •  (/ii  —  hf)  (23) 

Oxidant  pre-heating  :  0OX  =  Gox  •  (/12  —  ^4)  (24) 

Heat  recovered  :  0th  =  1^8  —  ^7 1  (25) 

BoP  thermal  balance  :  Wei  =  /*i  +  /i2  —  ^8  —  0th  —  0d_i  — 

0d 2  —  0d shell  +  ^yNif—AH)  (26) 


•  model  does  not  account  for  direct  oxidation  of  CO  at  cell 
anode; 

•  gas  temperatures  at  outlet  conditions  from  the  reformer 
and  the  stack  are,  respectively,  reforming  temperature  and 
operating  cell  temperature; 

•  internal  water  and  thermal  management. 


The  model  solves  a  non-linear  system  of  equations.  The  com¬ 
position  of  the  anodic  recycling  gas  is  then  calculated: 


< 


CO/  +  x  -  y  \  (H2i  +  3x  +  y  —  Z 


Uj  +  2x 


ni  +  2x 


i/n  V  'O  ~ r~  __  '  V _  1  ' 

Reforming  /  CH4/  _  x  \  /  H20/  -  X  -  y  +  Z 


^-Pshift 


Hi  +  2x 

'  H2/  +  3x  +  y  -  z 

Hi  +  2x 
CO  j  +  v  —  y' 


Hi  +  2x 
CO  2i  +  y 
Hi  +  2x 


Hi  +  2x 

z  =  Uf  ■  (3x  +  y  +  H2) 


H20/  -  x  -  y  +  z 
Hi  +  2x 


Pstack 


(14) 


Evaluated  x,  y  and  z  it  is  possible  to  define  the  anode  exhaust 
composition.  Most  of  the  considered  equations  are  listed  in 
Table  2. 

All  the  described  chemical  and  electric  calculations  need 
the  knowledge  of  the  stack  temperature.  A  tentative  stack 
temperature  is  assumed  and  the  cell  characteristics  evaluated. 
Then,  the  enthalpy  variation  inside  the  stack  is  calculated  and 
a  new  stack  temperature  obtained.  The  process  is  iterated  until 
the  satisfaction  of  a  prescribed  tolerance. 

The  results  obtained  with  the  model  have  been  compared 
with  the  experimental  data  of  the  CHP  100  acquired  during  its 
first  installation  at  EDB/ELSAM  (Westervoort),  data  which 
we  have  used  for  the  first  validation  of  the  model. 


5.  Methodology  and  computer  experimental  design 

Several  variables  can  be  used  to  characterize  the  stack 
operation,  and  they  can  be  considered  as  dependent  variables. 
These  dependent  variables  of  the  stack  are  influenced  by 


several  other  independent  variables  (in  the  following  defined 
factors ),  and  to  analyse  their  influence  an  experimental 
design  methodology  is  fundamental.  In  this  work,  the  exper¬ 
imental  tests  have  been  performed  in  the  form  of  a  computer 
experimental  session,  based  on  the  simulation  model  of  the 
stack;  the  measurement  uncertainties  are  simulated  with  per¬ 
turbation  imposed  to  the  independent  variables.  Therefore, 
a  full  factorial  design  of  experiments  has  been  employed:  an 
output  (dependent)  variable  could  be  correlated  to  the  input 
(independent)  variables  (or  factors)  pointing  out  its  effect. 
So,  it  is  possible  to  find  which  factors  deserve  a  particular 
attention:  for  example,  in  the  case  of  a  significant  factor  a 
further  analysis  could  be  done  or  more  accuracy  could  be 
necessary;  on  the  contrary,  the  factors  having  a  negligible 
effect  could  even  be  not  considered  in  subsequent  analysis. 
Also,  the  combined  effects  among  the  factors  can  be  assessed. 

Once  the  k  input  factors  have  been  chosen,  it  is  funda¬ 
mental  to  determine  their  range  of  variation  (the  experimental 
domain).  If  the  experiments  are  performed  only  at  the  extreme 
of  the  range  of  variation  of  the  factor  (that  is,  only  two  levels 
of  the  factor  are  considered),  2k  experiments  are  required,  and 
each  different  combination  of  different  factor  levels  is  called 
treatment.  If  the  treatment  is  replicated  then  an  analysis  of 
variance  can  be  performed. 

Concerning  the  adopted  symbolism,  each  factor  is  labelled 
with  a  capital  letter;  to  describe  the  treatments,  a  lowercase 
letter  is  used  when  the  factor  is  at  the  upper  level,  while  the 
letter  is  omitted  when  it  is  at  its  lower  level;  when  all  the 
factors  are  at  the  lower  level,  the  notation  (L)  is  used.  To 
analyse  the  effects  of  the  main  factors  the  well-known  Yates’ 
Technique  has  been  employed:  the  treatment  combinations 
must  be  therefore  written  down  in  standard  form  [48,49]. 
The  effects  of  the  factors  (indicated  with  the  letter  w),  are 
distinguished  in  main  effects  (wa,  wb,  . . .),  and  interaction 
effects  (wab,  wAc,  •  •  •)• 

In  the  experiment  session  discussed  in  the  paper,  five  fac¬ 
tors  have  been  chosen,  with  three  test  replications.  Therefore, 
a  total  of  25  x  3  x  2  =  192  computer  experiments  have  been 
necessary  to  obtain  all  the  possible  effects  and  interactions. 
The  five  factors  chosen  have  been:  air  utilization  factor  Uox, 
fuel  utilization  factor  t/p,  fuel  pre-heating  temperature  7p,  air 
pre-heating  temperature  Tox,  fraction  of  anode  gas  recycling 

P- 

The  range  of  the  factors  is  shown  in  Table  3. 

In  order  to  apply  the  Yates’  method  and  to  evaluate  the 
main  and  the  interaction  effects  of  the  factors,  we  need 
repetitions  on  the  experiments.  Therefore,  randomized 

Table  3 


Factors  analysed  (each  factor  is  labelled  with  a  capital  letter)  and  their  levels 


Factor 

Description 

Fower  level  (— ) 

Upper  level  (+) 

A 

Gox 

0.2 

0.3 

B 

Gf 

0.7 

0.9 

C 

Tv  (°C) 

200 

400 

D 

T)X  (°C) 

700 

850 

E 

P 

0.6 

0.7 
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Table  4 


Randomized  uncertainties  used  to  modify  the  factor  values 


Factor 

Uncertainty 

C/ox  (%) 

5 

UF  (%) 

2 

Tv  (°C) 

±1 

Fox  (°C) 

±1 

P(%) 

5 

uncertainties  on  the  selected  factors  (representing  indepen¬ 
dent  variables  of  the  model,  and  set  points  for  the  real  opera¬ 
tion)  and  on  the  current  (independent  variable  of  the  model, 
and  set  point  for  the  real  operation)  have  been  introduced.  The 
aim  was  to  simulate  a  real  experimental  measurement,  and 
the  randomized  uncertainties  take  into  account  the  errors  of  a 
real  instrument.  The  evaluated  effects  are  related  to  the  signif¬ 
icance  of  a  factor  macroscopic  variation  on  some  dependent 
variables,  and  not  to  the  effect  of  a  factor  measurement  error 
on  the  dependent  variables.  The  randomized  uncertainties 
used  to  modify  the  factors  values  are  shown  in  Table  4. 

Six  outputs  are  investigated  in  this  work:  the  electric  Wei 
and  recovered  heat  power  d>th  outputs;  the  fuel  mass  flow 
Gp;  the  single  cell  voltage  Vc;  the  stack  temperature  rstack; 
the  steam- to-carbon  ratio  S/C. 

To  reduce  the  number  of  experiments,  only  two  operat¬ 
ing  conditions  have  been  analysed  at  two  different  current 
density  values  (0.2  and  0.4  A  cm-2),  introducing  also  a  ran¬ 
domized  uncertainty  on  the  current  data  (1%,  simulating  the 
uncertainty  in  the  current  acquisition).  The  factorial  analy¬ 
sis  was  done  at  the  indicated  two  current  densities,  for  every 
dependent  variable.  These  values  of  current  density  have  been 
chosen  because  higher  values  are  not  realistic  for  the  analysed 
stack.  Actually,  even  the  value  of  0.4  A  cm-2  is  very  high, 
but  some  reports  states  that  it  could  be  reached  at  operating 
temperatures  higher  than  900  °C  with  a  cell  voltage  around 
0.47  V  [43]. 

The  computer  experimental  campaign  has  been  carried  out 
considering  the  EOS-CHP100  fed  with  dry  methane. 


6.  Results  and  discussion 

6.1.  Electric  power 

6.1.1.  Considerations  about  main  and  interaction  effects 
The  main  and  interaction  effects  on  electric  power  are 

shown  in  Fig.  3  (low  current  density)  and  Fig.  4  (high  current 
density). 

Factors  with  no  significant  effects:  An  interesting  notation 
is  that  the  fuel  preheating  temperature  has  no  effect  on  the 
electric  power  produced;  therefore,  the  fuel  preheating  could 
be  avoided.  Moreover,  even  the  anode  gas  recycling  fraction 
shows  no  significant  effects,  and  therefore  it  could  be  fixed 
at  the  value  indicated  by  the  plant  manuals. 


Uox  at  high  level:  Positive  main  effect,  because  of  better 
performances  in  terms  of  internal  thermal  management  and 
single  cell  electric  efficiency,  due  to  increase  of  cell  temper¬ 
ature  and  decrease  of  cell  overpotentials. 

Up  at  high  level:  Negative  main  effect,  due  to  decrease  of 
hydrogen  partial  pressure  at  cell  anode,  decrease  in  Nernst 
voltage. 

Tox  at  high  level:  Positive  main  effect,  due  to  increase  of 
cell  temperature  and  decrease  of  cell  overpotentials. 

Up  at  high  level,  Tox  at  high  level:  Negative  effect,  because 
of  prevailing  negative  effect  of  Gp  at  high  level. 

Up  at  high  level,  f  at  high  level:  Treatment  combinations 
with  high  fuel  utilization  and  high  anode  gas  recycling  frac¬ 
tion  have  a  minor  negative  effect,  because  of  the  negative 
effect  of  Gp  at  high  level  and  of  the  excess  water  in  case  of 
higher  fraction  of  anode  gas  recycling. 

Main  and  interaction  effects  at  high  cell  current:  At  high 
cell  currents  the  main  effects  of  the  factors  are  increased 
because  of  bigger  mass  flows  inside  the  stack;  moreover, 
at  high  cell  current  the  fraction  of  anode  gas  recycling  has 
a  main  minor  negative  effect  on  electric  power,  because  of 
the  excess  water  concentration  at  the  anode  side  amplified 
at  high  current. 

6.1.2.  Analysis  of  the  effects  of  the  significant  factors 

The  effect  of  the  air  utilization  factor  is  shown  in  Fig.  5. 

The  curves  with  high  air  utilization  factor  show  higher 
values  of  stack  electric  power  at  various  current  density  val¬ 
ues.  A  higher  Gox  means  a  lower  air  mass  flow,  and  therefore 
a  higher  stack  equilibrium  temperature  (at  fixed  conditions). 
The  higher  stack  temperature  causes  an  increase  of  the  cell 
voltage:  it  is  possible  to  obtain  especially  a  reduction  of  the 
activation  and  ohmic  overvoltages. 

The  effects  are  increased  at  high  current  densities,  because 
the  air  mass  flows  are  greater  and  this  amplifies  the  effects  of 
changing  the  air  utilization  factor. 

A  factor  causing  significant  effects  on  the  stack  electric 
power  is  the  air  pre-heating  temperature  Tox ,  as  shown  in 
Fig.  6. 

The  increase  of  the  air  pre-heating  temperature  causes 
an  increase  of  the  stack  equilibrium  temperature,  and  con¬ 
sequently  an  increase  of  the  cell  voltage:  it  is  possible  to 
obtain  especially  a  reduction  of  the  activation  and  ohmic  over¬ 
voltages.  The  effects  are  increased  at  high  current  densities, 
because  the  air  mass  flows  are  greater  and  this  amplifies  the 
effects  of  changing  the  air  utilization  factor. 

Another  factor  causing  significant  effects  on  the  stack 
electric  power  is  the  fuel  utilization  factor  Gp,  as  shown  in 
Fig.  7. 

An  increase  of  the  factor  Gp  determines  a  reduction  of 
the  stack  electric  power;  this  is  due  to  a  reduction  of  the  cell 
voltage,  caused  by  an  increase  of  the  molar  fraction  of  H2O  in 
the  anode  mixture,  with  a  consequent  reduction  of  the  molar 
fraction  of  H2  on  the  electrode  surface,  as  shown  in  Table  5. 

As  usual,  the  effects  are  emphasized  at  higher  currents  due 
to  the  higher  mass  flows. 
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Fig.  3.  Computed /-value  and  effects  of  the  parameters  on  stack  electric  power  and  heat  recovered  at  i  =  0.2  A  cm 


-2 


M.  Cali  et  al.  /  Journal  of  Power  Sources  156  (2006)  400-413  407 


Computed  f  on  Electric  Power 


ic=0.2  A/cm2 - alfa=4.98 


®  ^  ®  ®&  6  £  e  <?  ■JjP  o®^o%o®  o®  ^VVV^VVV® 


O'  v  v  v  <2r 

Treatment  Combinations 


Factors's  Effect  on  Electric  Power:  ic=0.4  A/cm‘ 


10 

9 

8 

7 

6 

-  5 
4 
3 
2 
1 
0 


Treatment  Combinations 

Factors's  Effect  on  Heat  Recovered:  ic=0.4  A/cm2 


10 

8 

6 

4 

2 

co 

o 

ig  0 

LU 

-2 

-4 

-6 

-8 

-10 


r- 

— 

n 

0 

b  n  FI  Cl  n 

n 

(0  t 

:T1JT¥T! 

-t - 1 - 1 - 1 - 1 - rrj-r  n  r 

T5  TJ  *0  TJ  D  ^ 

”  ■“  1  °  8  J  3 

CO 

i  r.  - 1 - 1 - 1 - r  i— i  i - 1 - 1  i - rr  1 1 - 1 - 1 - r«-< 

iQ<D<Da)Q)<D<u  a>  a>  M  a>  <d  a  a> 

a  JD-Q  O  O  O  o  -o  T3  -D  ■0  T3  13  C  "O 

CO  CO  -Q  -Q  co  -Q  “  o  o  c_  o 

LI  co  a)  co  (D  lu  ^ 

Treatment  Combinations 


Computed  f  on  Heat  Recovered 


ic=0.2  A/cm2 - alfa=4.98 


Fig.  4.  Computed /-value  and  effects  of  the  parameters  on  stack  electric  power  and  heat  recovered  at  i  =  0.4  A  cm 


408  M.  Cali  et  al.  /  Journal  of  Power  Sources  156  (2006)  400-413 


M.  Cali  et  al.  /  Journal  of  Power  Sources  156  (2006)  400-413 


409 


Effect  of  Uox  on  Electric  Power  (DC)  at  various  current  density  values 
Tf=25°C;  TOx=850°C;  Uf=0.85;  ps,ack=1bar 


~+~  Uox=0.2  Uox=0.25  *  Uox=0.3 


i  (A/cm2) 


Fig.  5.  Effect  of  the  air  utilization  factor  on  stack  electric  power. 


Effect  of  Tox  on  Electric  Power  (DC)  at  various  current  density  values 
Tf=25°C;  Uf=0.85;  Uox=0.3;  pslack=1bar 


Tox=700°C  Tox=750°C  Tox=800°C  -*-Tox=850°C 


i  (A/cm2) 


Fig.  6.  Effect  of  the  air  pre-heating  temperature  on  stack  electric  power. 


Effect  of  UF  on  Electric  Power  (DC)  at  various  current  density  values 
Tf=25°C;  TOx=850°C;  Pstack=1  bar;  UOX=0.25 


UF=0.85  UF=0.9 


i  (A/cm2) 


Fig.  7.  Effect  of  the  fuel  utilization  factor  on  stack  electric  power. 


Table  5 

Effect  of  Up  on  the  molar  fraction  of  H2  on  the  anode  mixture  and  consequent 


effect  on  the  cell  open  circuit  voltage 


Treatment  combination 

^h2o 

^h2 

Jco 

yco2 

Erev 

(L) 

0.4858 

0.148 

0.06 

0.3062 

0.864 

b 

0.6114 

0.0452 

0.025 

0.3184 

0.777 

Combined  effect  of  Tox  and  Uox  on  Electric  Power  (DC)  at  various  current 

density  values 


Uox=0.2  Uox=0.25  *  Uox=0.3 


Tox  (°C) 


Fig.  8.  Combined  effect  of  the  air  utilization  factor  and  the  air  pre-heating 
temperature  on  stack  electric  power. 

It  is  interesting  to  analyze  the  effects  of  some  factor  com¬ 
binations.  The  factor  combinations  with  a  significant  effect 
on  the  stack  electric  power  are  Uox-Tox  and  U^-Tox.  The 
effect  of  the  combination  of  Uox-Tox  on  the  stack  electric 
power  is  shown  in  Fig.  8. 

We  observe  that  the  effects  of  the  combination  are  impor¬ 
tant  especially  at  high  current  values.  In  fact,  the  effects  of 
the  air  utilization  factor  are  emphasized  due  to  the  higher 
mass  flows.  Moreover,  the  effect  of  the  air  utilization  factor 
is  important  if  the  air  pre-heating  temperature  is  low,  while 
at  high  values  of  the  air  pre-heating  temperature  the  effect  of 
the  air  utilization  factor  is  hardly  reduced.  In  fact,  the  positive 
effect  of  a  high  air  utilization  factor  on  the  stack  equilib¬ 
rium  temperature  is  evident  if  the  incoming  air  is  slightly 
pre-heated,  while  if  the  incoming  air  has  been  already  pre¬ 
heated  at  high  temperature  the  effect  of  the  air  utilization 
factor  on  the  stack  temperature  is  reduced. 

Finally,  the  effect  of  the  combination  of  U^-Tox  on  the 
stack  electric  power  is  shown  in  Fig.  9. 

Considering  a  combination  of  high  values  of  I/p  and  Tox, 
the  reduction  of  the  stack  electric  power  due  to  a  high  fuel  uti¬ 
lization  factor  could  be  compensated  by  a  higher  pre-heating 


Combined  effect  of  Tox  and  U  F  on  Electric  Power  (DC)  at  various  current 

density  values 

-^UF=0.85  -*-UF=0.9 


Tax  (°  C) 


Fig.  9.  Combined  effect  of  the  fuel  utilization  factor  and  the  air  pre-heating 
temperature  on  stack  electric  power. 
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of  the  incoming  air.  This  effect  is  particularly  emphasized  at 
high  current  densities,  where  the  curve  slope  is  accentuated. 

6.2.  Recovered  heat  power 

6.2.1.  Considerations  about  main  and  interaction  effects 
The  main  and  interaction  effects  on  heat  recovered  are 

shown  in  Fig.  3  (low  current  density)  and  Fig.  4  (high  current 
density). 

Factors  with  no  significant  effects:  An  interesting  notation 
is  that  the  fuel  preheating  temperature  has  no  effect  on  the 
recovered  heat  power;  therefore,  the  fuel  preheating  could 
be  avoided. 

U ox  at  high  level:  Positive  main  effect,  due  to  increase  of 
exhaust  gas  temperature  and  decrease  of  heat  used  for  oxi¬ 
dant  preheating. 

Uf  at  high  level:  Negative  main  effect,  due  to  decrease  of 
heat  recovered  from  the  depleted  fuel  in  the  post-combustion 
chamber. 

Tox  at  high  level:  Negative  main  effect,  due  to  better  electric 
performances  of  the  stack  (and  reduction  of  heat  produced 
by  irreversibility)  and  increase  of  oxidant  preheating. 
ft  at  high  level:  Negative  effect,  an  increase  of  ft  reduces  the 
cell  efficiency  (and  therefore  increases  the  heat  generated 
by  irreversibility)  but  at  the  same  time  the  fuel  fraction  at 
the  post-combustor  is  reduced,  and  this  is  the  prevailing 
factor  explaining  the  negative  effect;  this  negative  effect  is 
increased  with  a  high  fuel  utilization  factor  (combination 
be). 

Main  and  interaction  effects  at  high  cell  current:  At  high 
cell  current  the  main  effects  of  the  factors  are  emphasized. 

6.2.2.  Analysis  of  the  effects  of  the  significant  factors 

The  effect  of  the  air  utilization  factor  is  shown  in  Fig.  10. 
With  an  increase  of  the  air  utilization  factor,  the  heat  recov¬ 
ered  for  cogeneration  increases.  This  is  due  to  several  fac¬ 
tors:  higher  stack  equilibrium  temperature,  and  consequently 
higher  temperature  of  the  stack  exhaust  flow;  moreover,  the 
heat  flow  used  to  preheat  the  oxidant  is  reduced  because  the 

Effect  of  Uox  on  Heat  Recovered  at  various  current  density  values 
Tf=25°C;  T ox=850°C;  Uf=0.85;  pstack="lbar 


^  Uox=0.2  Uox=0.25  Uox=0.3 


i  (A/cm2) 


Fig.  10.  Effect  of  the  air  utilization  factor  on  heat  recovered. 


Effect  of  Tox  on  Heat  Recovered  at  various  current  density  values 
Tf=25°C;  Uf=0.85;  Uox=0.3;  ps,ack=1bar 

“•~Tox=700°C  Tox=750°C  *  Tox=800°C  -^Tox=850°C 


Fig.  11.  Effect  of  the  air  pre-heating  temperature  on  heat  recovered. 

inlet  cathode  mass  flow  is  lower  at  lower  air  utilization  factor. 
Even  in  this  case,  the  effects  are  increased  at  high  current  den¬ 
sities,  because  the  air  mass  flows  are  greater  and  this  amplifies 
the  effects  of  changing  the  air  utilization  factor. 

A  factor  causing  significant  effects  on  the  heat  recovered 
is  the  air  pre-heating  temperature  Tox ,  as  shown  in  Fig.  11. 

The  increase  of  the  air  pre-heating  temperature  has  a  neg¬ 
ative  effect  on  the  heat  recovered  for  cogeneration.  This  is 
due  to  the  fact  that  the  heat  used  to  preheat  the  incoming  air 
increases  with  the  Tox  factor,  and  therefore  the  heat  recov¬ 
ered  for  cogeneration  is  reduced.  This  happens  even  if  the 
Tox  increase  has  a  positive  effect  on  the  stack  temperature 
and  the  exhaust  temperature  exiting  the  stack.  As  usual,  the 
effects  are  emphasized  at  higher  currents  due  to  the  higher 
mass  flows. 

The  effect  of  the  fuel  utilization  factor  t/p  on  the  heat 
recovered  is  shown  in  Fig.  12. 

It  is  interesting  to  observe  that  the  heat  recovered  for 
cogeneration  increases  at  low  fuel  utilization  factor  values.  In 
fact,  with  high  f/p  the  fuel  mass  flow  at  the  post-combustor  is 
reduced,  reducing  the  heat  generated  by  post-combustion, 
and  therefore  the  stack  exhaust  temperature  is  lower.  To 
extract  a  higher  heat  flow  from  the  stack  exhaust  is  there¬ 
fore  better  to  operate  with  lower  fuel  utilization  factor  values. 

Effect  of  Uf  on  Heat  Recovered  at  various  current  density  values 
Tf=25°C;  TOx=850°C;  pstack=1  bar;  U  ox=0-25 


UF=0.8  5  UF=0.9 


i  (A/cm2) 

Fig.  12.  Effect  of  the  fuel  utilization  factor  on  heat  recovered. 
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Combined  effect  of  Toxand  Uoxon  Heat  Recovered  at  various  current  density 

values 


-*~Uox=0.2  -*-Uox=0.25  Uox=0.3 


Fig.  13.  Combined  effect  of  the  air  utilization  factor  and  the  air  pre-heating 
temperature  on  heat  recovered. 

This  happens  even  if  at  high  I/p  the  heat  generated  by  irre¬ 
versibility  increases  (as  we  noted,  in  this  case  the  cell  voltage 
is  reduced),  determining  a  higher  stack  equilibrium  tempera¬ 
ture,  but  this  is  not  sufficient  to  assure  an  increase  of  the  heat 
recoverable  from  the  exhaust  flow. 

It  is  interesting  to  analyze  the  effects  of  some  factor  combi¬ 
nations.  The  factor  combinations  with  a  significant  effect  on 
the  heat  recovered  for  cogeneration  are  Uox-Tox  and  U^-Tox. 
The  effect  of  the  combination  of  Uox-Tox  on  the  heat  recov¬ 
ered  is  shown  in  Fig.  13. 

We  observe  that  the  combination  determines  different 
effects  at  low  or  high  current  densities. 

At  fixed  Tox,  a  higher  air  utilization  factor  causes  an 
increase  of  the  stack  equilibrium  temperature,  with  a  con¬ 
sequent  increase  of  the  cell  voltage  and  reduction  of  the 
heat  generated  by  irreversibility;  at  the  same  time,  there  is 
a  reduced  mass  flow  of  incoming  air  which  has  to  be  pre¬ 
heated  using  the  heat  flow  taken  from  the  stack  exhaust. 

At  low  current  densities,  a  high  air  utilization  factor  deter¬ 
mines  a  higher  value  of  stack  electric  power  for  every  value  of 
air  pre-heating  temperature.  This  could  be  due  to  the  reason 
that  in  these  operating  conditions  (low  current,  and  therefore 
low  stack  irreversibility)  the  prevailing  effect  is  the  reduc¬ 
tion  of  the  incoming  air,  which  has  to  be  pre-heated,  with  a 
consequent  increase  of  the  heat  recoverable  for  cogeneration. 

At  high  current  densities,  the  behavior  is  more  complex. 
The  operation  irreversibility  (and  therefore  the  heat  gener¬ 
ated  by  the  stack)  is  higher  than  at  low  current  densities,  and 
the  role  of  the  irreversibility  generated  heat  is  emphasized. 
Moreover,  as  we  have  observed  in  Fig.  12,  the  stack 
equilibrium  temperature  variation  due  to  the  variation  of 
the  air  utilization  factor  is  emphasized  at  low  values  of  air 
pre-heating  temperature.  Therefore,  if  a  low  Tox  is  fixed,  at 
high  currents  is  convenient  to  send  to  the  stack  a  greater  air 
mass  flow  (low  air  utilization  factor)  in  order  to  exploit  the 
higher  irreversibility  generated  heat  to  obtain  an  increase  of 
the  heat  recovered  for  cogeneration.  Conversely,  if  a  high 
Tox  is  fixed,  the  behavior  is  the  same  observed  at  low  current 
densities:  in  order  to  increase  the  heat  recovered  for  cogen¬ 


Combined  effect  of  Tox  and  UF  on  Heat  Recovered  at  various  current  density 

values 


UF=0.8  5 UF=0.9 


T0x(°C) 


Fig.  14.  Combined  effect  of  the  fuel  utilization  factor  and  the  air  pre-heating 
temperature  on  heat  recovered. 

eration  it  is  convenient  to  operate  with  a  higher  air  utilization 
factor. 

Finally,  the  effect  of  the  combination  of  Up-Tox  on  the 
heat  recovered  is  shown  in  Fig.  14. 

Considering  a  combination  of  high  values  of  f/p  and  Tox, 
the  effect  is  the  reduction  of  the  heat  recovered  for  cogen¬ 
eration.  The  negative  effect  is  emphasized  at  high  current 
densities. 


7.  Conclusions 

In  this  paper,  the  computer  experimental  analysis,  through 
a  statistical  methodology  (ANOVA  with  Yate’s  method),  of 
the  single  and  combined  effects  of  the  main  stack  operation 
independent  variables  on  the  cogenerative  (electricity  and 
heat)  performance  of  the  EOS-CHP  100  SOFC  Field  Unit 
(Siemens)  with  methane  feeding  has  been  developed.  The 
results  can  be  summarized  as  follows: 

•  The  fuel  preheating  temperature  has  no  effect  on  both  the 
electric  power  produced  and  the  heat  recovered;  therefore, 
the  fuel  preheating  could  be  avoided. 

•  The  anode  gas  recycling  fraction  /3  shows  no  significant 
effects  on  the  electric  power,  while  it  shows  a  negative 
effect  on  the  heat  recovered,  because  the  fuel  fraction  at 
the  post-combustor  is  reduced;  therefore  it  could  be  fixed 
at  the  value  indicated  by  the  plant  manuals. 

•  U ox  at  high  level :  It  shows  positive  main  effect  on  the 
electric  power,  because  of  better  performances  in  terms  of 
internal  thermal  management  and  single  cell  electric  effi¬ 
ciency,  due  to  increase  of  cell  temperature  and  decrease  of 
cell  overpotentials;  at  the  same  time,  it  shows  positive  main 
effect  on  the  heat  recovered,  due  to  increase  of  exhaust  gas 
temperature  and  decrease  of  heat  used  for  oxidant  preheat¬ 
ing. 

•  Up  at  high  level :  It  shows  negative  main  effect  on  the  elec¬ 
tric  power,  due  to  decrease  of  hydrogen  partial  pressure  at 
cell  anode  and  consequent  decrease  in  Nernst  voltage;  at 
the  same  time,  it  shows  a  negative  main  effect  on  the  heat 
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recovered,  due  to  decrease  of  heat  recovered  from  wasted 
fuel  in  post-combustion  chamber. 

•  Tox  at  high  level :  It  shows  positive  main  effect  on  the  elec¬ 
tric  power,  due  to  increase  of  cell  temperature  and  decrease 
of  cell  overpotentials;  at  the  same  time,  it  shows  a  negative 
main  effect  on  the  heat  recovered,  due  to  better  electric 
performances  of  the  stack  (and  reduction  of  heat  produced 
by  irreversibility)  and  increase  of  oxidant  preheating. 

•  At  high  cell  current  densities  the  main  effects  of  the  factors 
are  emphasized  because  of  the  bigger  mass  flows  inside  the 
stack. 
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